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LONG  TERM  GOALS 

The  long  terms  goals  of  the  project  are  to  develop  analytical  or  efficient  numerical  schemes  for 
determining  the  RE  signal  received  over  a  rough  ocean  surface  and  in  atmospheric  ducting  conditions. 

OBJECTIVES 

The  objectives  of  the  proposed  work  are  to  develop  techniques  for  treating  radiowave  propagation  over 
an  electrically  large  rough  surface.  Both  analytical  and  numerical  techniques  were  looked  at  in  this 
project. 

APPROACH 

One  graduate  Ph.D.  student,  Zhiguo  Eai,  was  fully  supported  through  the  grant  to  help  carry  out  the 
current  research  work.  He  is  expected  to  complete  his  PhD  by  January  2007. 

Significant  advances  have  been  made  in  understanding  the  basic  phenomena  by  means  of  experimental 
investigations  of  electromagnetic  scattering  from  the  ocean  surface  [1,2]  and  in  controlled  wave  tank 
experiments  [3,4].  Meanwhile  a  large  number  of  publications  on  the  theoretical  treatment  of  the 
general  problem  are  available  in  literature  [5,6].  However  the  theory  is  still  far  from  completely 
developed  as  almost  all  the  existing  scattering  models  are  only  valid  for  small  roughness  and  most  of 
them  are  expected  to  fail  at  low  grazing  incidence.  Rigorous  numerical  methods  have  been  widely  used 
over  the  past  several  decades  to  overcome  the  limitations  of  analytical  theories  and  to  provide 
validation  for  approximate  scattering  models  [7,  8].  Eor  one-dimensional  problems,  available 
computational  resources  now  allow  non-grazing  problems  to  be  readily  solved  using  any  of  a  number 
of  standard,  well  established  approaches.  The  major  remaining  challenge  in  the  development  of 
numerical  methods  is  the  low  grazing  angle  scattering  problem.  For  example  the  standard  method  of 
moment  has  not  proven  well  suited  for  applications  at  low  grazing  angles  as  surfaces  of  increasingly 
large  size  are  required  for  simulations  when  the  direction  of  incidence  approaches  horizontal.  Larger 
surfaces  mean,  in  turn,  larger  matrix  sizes  and  higher  computational  and  storage  expenses.  Furthermore 
surfaces  of  finite  length  are  used  in  numerical  simulations  leading  to  nonphysical  edges  that  can 
introduce  strong  diffraction  in  the  scattered  field. 

In  our  current  work,  a  formulation  based  on  full  wave  integral  equation  was  developed.  Our  primary 
goal  is  to  study  the  statistical  properties  of  the  scattering  amplitude  and  the  scattered  fields  and  develop 
a  rough  surface  scattering  model  that  is  especially  suitable  for  scattering  from  a  very  rough  surface  at 
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low  grazing  angles.  The  rough  surfaee  is  assumed  to  be  periodic  over  lengths  much  larger  than  the 
correlation  length  to  avoid  artificial  edge  effects.  An  ensemble  of  periodic  surfaces  with  appropriate 
statistical  properties  was  considered  to  generate  the  statistical  properties  of  the  field.  The  period  of  the 
surface  is  chosen  to  be  large  enough  so  that  there  are  many  peaks  and  valleys  of  the  height  profile 
within  one  period  and  the  results  obtained  from  periodic  surfaces  are  a  good  predictor  for  the  infinite- 
surface  case  of  coherent/incoherent  power,  etc.  We  intend  to  answer  questions  such  as  ‘whether  the 
coherent  reflection  coefficient  is  adequate  to  describe  the  scattering  phenomena’  and  ‘whether  the 
coherent  component  continues  to  dominate  for  large  roughness  and  low  grazing  angles  and  at  what 
point  the  diffuse  component  takes  over’  as  they  have  not  been  addressed  in  the  literature  definitively. 

To  start  we  consider  two-dimensional  problems  with  Dirichlet  boundary  conditions  where  a  one¬ 
dimensional  periodic  rough  surface  is  illuminated  by  a  down-propagating  plane  wave  and  the  total 
electric  field  on  the  rough  surface  vanishes.  The  starting  point  in  our  formulation  is  the  Helmholtz 
formula  leading  to  a  surface  integral  equation  for  the  induced  surface  current.  The  periodic  boundary 
condition  is  then  analytically  incorporated  resulting  in  a  surface  integral  over  only  one  single  period 
instead  of  over  an  infinite  domain.  The  integral  equation  is  then  discretized  through  the  method  of 
moments  and  solved  using  the  generalized  conjugate  residue  algorithm  with  the  help  of  the  fast 
multipole  method  [9].  Significant  progress  has  already  been  made  in  the  development  of  this 
formulation  and  some  preliminary  results  have  been  obtained. 

Usually  it  is  convenient  to  express  the  received  signal  in  terms  of  the  propagation  factor,  defined  as  the 
ratio  of  the  total  field  to  the  incident  field  due  to  the  direct  path: 

F  F  F  F 
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The  propagation  factor  is  further  separated  into  coherent  and  incoherent  components.  The  coherent 
component  is  the  power  of  average  field  and  is  expressed  as  =  \\Fj\  .  The  incoherent  component  is 

the  mean  of  the  variance  of  E  and  is  expressed  as  In  these  definitions,  denotes 

ensemble  average  over  many  realizations. 

All  numerical  simulation  results  below  were  performed  at  a  frequency  of  3  GHz  (wavelength  /I  =  0.1 
m).  The  plane  wave  was  incident  at  a  grazing  angle  of  5.0  degrees.  Electric  field  was  computed  on  a 
vertical  line.  Rough  surfaces  were  generated  using  a  zero-mean  Gaussian  height  statistics  and  Pierson- 
Moskowitz  spectrum  which  is  a  function  of  the  wind  speed  U  at  a  height  of  19.5  m.  The  step-size 
along  x-direction  is  /1/8  and  the  overall  surface  length  is  2000/1  (16,000  unknowns).  Ensemble  averages 
were  taken  over  1000  realizations. 

Eigures  1  to  4  show  the  coherent/incoherent  power  for  U=2,5,  10,  15  m/s,  respectively.  It  is  seen  that 
the  difference  between  the  coherent  and  incoherent  parts  decreases  as  the  surface  becomes  rougher. 
Eigures  5  to  8  show  the  histograms  of  the  propagation  factor  (real  part,  imaginary  part  and  magnitude) 
superimposed  by  analytical  Gaussian  PDEs  at  a  fixed  height  (z  =  5  m)  for  U  =2,5,  10,  15  m/s, 
respectively.  The  scattered  field  closely  follows  a  complex  Gaussian  random  process  with  independent 
real  and  imaginary  parts.  The  correlation  coefficients  between  the  real  part  and  imaginary  part  are 
0.036,  0.063,  0.001,  and  0.082,  respectively. 


Figure  1.  Coherent  and  incoherent  power  for  U  =  2  m/s. 
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Figure  2.  Coherent  and  incoherent  power  for  U=  5  m/s. 
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Figure  3.  Coherent  and  incoherent  power  for  U  =  10  m/s. 
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Figure  4.  Coherent  and  incoherent  power  for  U  =  15  m/s. 


Figure  5.  Histograms  (PDFs)  of  the  propagation  factor  at  z  =  5  m  for  U  =  2  m/s. 


Figure  6.  Histograms  (PDFs)  of  the  propagation  factor  at  z  =  5  m  for  U  =  5  m/s. 
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Figure  7.  Histograms  (PDFs)  of  the  propagation  factor  at  z  =  5  m  for  U  =  10  m/s, 
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Figure  8.  Histograms  (PDFs)  of  the  propagation  factor  at  z  =  5  m  for  U  =  15  m/s. 


IMPACT/APPLICATIONS 


The  results  presented  thus  far  and  successful  completion  of  the  ongoing  research  should  help  the  PE 
code  developers  in  generating  better  operational  codes  for  assessing  the  effects  of  ocean  roughness  on 
radar  detectability  and  surveillance. 
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